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ABSTRACT 

We present a new method to estimate the absolute ages of stellar systems. This method is based 
on the difference in magnitude between the main sequence turn-off (MSTO) and a well defined knee 
located along the lower main sequence (MSK). This feature is caused by the collisionally induced 
absorption of molecular hydrogen and it can be easily identified in near-infrared (NIR) and in optical- 
NIR color-magnitude diagrams of stellar systems. We took advantage of deep and accurate NIR images 
collected with the Multi-Conjugate Adaptive Optics Demonstrator temporarily available on the Very 
Large Telescope and of optical images collected with the Advanced Camera for Surveys Wide Field 
Camera on the Hubble Space Telescope and with ground-based telescopes to estimate the absolute age 
of the globular NGC 3201 using both the MSTO and the A(MSTO-MSK). We have adopted a new set 
of cluster isochrones and we found that the absolute ages based on the two methods agree to within 
one sigma. However, the errors of the ages based on the A(MSTO-MSK) method are potentially 
more than a factor of two smaller, since they are not affected by uncertainties in cluster distance 
or reddening. Current isochrones appear to predict slightly bluer («0.05 mag) NIR and optical-NIR 
colors than observed for magnitudes fainter than the MSK. 

Subject headings: globular clusters: individual (NGC3201) — stars: evolution — stars: fundamental 
parameters 
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1. INTRODUCTION 

The absolute ages of Galactic Globular Clusters 
(GGCs) is a crossroad of several astrophysical prob- 
lems (VandenBerg et al. 1996; Chaboyer 1998; Castellani 
1999). This parameter provides: i) a lower limit to the 
age of the universe (Buonanno et al. 1998; Stetson et al. 
1999; Gratton et al. 2003); ii) robust constraints on the 
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physics adopted in stellar evolutionary models (Salaris & 
Weiss 1998; Cassisi et al. 1999; VandenBerg et al. 2008; 
Dotter et al. 2008), and Hi) the chronology for the assem- 
bly of the halo, bulge and disk of the Milky Way (Carraro 
et al. 1999; Rosenberg et al. 1999; Zoccali et al. 2003; de 
Angeli et al. 2005). However, estimates of absolute GC 
ages are still hampered by uncertainties in the distance, 
extinction and the reddening law (Renzini 1991; Bono 
et al. 2008), the overall metallicity scale (Rutledge et al. 
1997; Kraft & Ivans 2003), as well as the detailed dis- 
tribution of heavy-element abundances (Asplund et al. 
2005), and photometric zero-points (Stetson 2000). 

Proposed ways of overcoming some of these thorny 
problems include the use either of a different clock (white 
dwarf cooling sequence, Richer et al. 2006), or different 
photometric systems (Stromgren bands, Grundahl et al. 
1998) or different diagnostics (luminosity function, Zoc- 
cali et al. 2000; Richer et al. 2008). However, these alter- 
native approaches to estimate the absolute ages of GCs 
can be impeded by limiting magnitude, by photometric 
accuracy, or by sample completeness. Simultaneous and 
self-consistent interpretation of optical and NIR photom- 
etry is another way to refine current distance and red- 
dening estimates. In fact, as shown in this study (see 
also Calamida et al. 2009) deep J, K photometry reveals 
the existence of a well-defined "knee" in the lower main 
sequence. The position of this feature is, at fixed chem- 
ical composition, essentially independent of the cluster 
age. Therefore, the difference either in magnitude or in 
color with the turn-off can provide robust estimates of 
the absolute age. 

Consider the GC NGC 3201, which has several inter- 
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csting properties. It is relatively nearby and it suffers 
from appreciable foreground reddening (/i = 13.32±0.06, 
E(B -V) = 0.30 ± 0.03, Piersimoni et al. 2002; fi = 
13.36 ± 0.06, E{B -V) = 0.25 ± 0.02, Layden & Saraje- 
dini 2003; Mazur et al. 2003). Accurate abundances esti- 
mates exist for both iron ([Fe/H]=-1.54 ± 0.10, Kraft & 
Ivans 2003; Covey et al. 2003) and a-elements ([a/Fe] — 
0.2-0.4, Pritzl et al. 2005). However, NGC 3201 is af- 
fected by field contamination and by differential redden- 
ing (Piersimoni et al. 2002; Kravtov et al. 2009). Owing 
to these problems, we still lack an accurate estimate of 
its absolute age (Layden & Sarajedini 2003). 

2. OBSERVATIONS AND DATA REDUCTION 

Near-Infrared (NIR) data (J, K s ) were collected in two 
observing runs with different pointings of the NIR cam- 
era SOFI (Field of View [FoV] ~ 5'x 5'; scale 0"29 /px) 
on the New Technology Telescope (NTT; ESO, La Silla). 
We obtained 132 J-band images with individual expo- 
sure times from 3 to 6 s and median seeing 0''8. We also 
collected 161 X s -band images with exposure times from 
10 to 12 s and median seeing 0"6. The total exposures 
per band are « 10 (J) and w 40 minutes (K s ). These 
images unevenly cover an area of « 20' x 18'across the 
cluster center (see Fig. 1). NIR data (J,K S ) were also 
collected with the NIR camera ISPI (FoV -10!5xl0^5; 
scale 0"30 /px) on the 4m Blanco telescope (CTIO; La 
Serena). We collected 15 J-band images, each of which 
was the co-added result of six individual 5 s exposures; 
these had median seeing 0"74 arcsec. We also obtained 
15 i^s-band images, consisting of eight co-added images 
with individual exposure times of 5 s, and median seeing 
0"75 arcsec. The total exposures per band are thus 450 s 
in J and 600s in K s . 

These data were supplemented with deep NIR data 
(J, K s ) from the Multi-Conjugate Adaptive Optics 
Demonstrator (MAD) temporarily available on the Very 
Large Telescope (VLT; ESO, Paranal). MAD is a proto- 
type instrument performing real-time correction for at- 
mospheric turbulence (Marchetti et al. 2006). Its in- 
frared camera employs a 2048x2048 Hawaii-II infrared 
detector with a scale of 0''028 /px for a total FoV of w 
1'. During the first Science Demonstration of MAD four 
l'xl' fields were observed in a region located in the S-W 
corner of NGC 3201. Five guide stars with visual mag- 
nitudes ranging from 11.8 to 12.9 distributed on a circle 
of 2 arcmin diameter concentric to the field were used. 
For each pointing we collected three J-band and five K s - 
band images of 240 sec (DIT^IO, NDIT=24). The seeing 
during the observations changed from 0.6 to 0.9 arcsec 
(J-band) and from 0.8 to 1.3 arcsec (if s -band). The full 
width at half-maximum (FWHM) measured on the im- 
ages ranges from 0.07 to 0.10 arcsec. We ended up with a 
catalog containing —54,100 stars with at least one mea- 
surement in each of two different bands. The top panel 
of Fig. 1 shows the spatial distribution of the different 
NIR datasets 17 . 

The V and / data used in this investigation come from 
the database of original and archival observations that 
have been collected, reduced, and calibrated by one of 
the authors (PBS) in an ongoing effort to provide homo- 

17 Following Harris (1996) we assumed for the center of the clus- 
ter the coordinates: a = 10 h 17 m 36"8, <5=-46°24'40" 



geneous photometry on the Landolt (1992) photometric 
system for a significant fraction of GGCs (Stetson 2000). 
For NGC 3201, we currently have a catalog consisting of 
—180,000 stars with at least two measurements in each 
of B, V, and /; these span an area with extreme dimen- 
sions on the sky of ~40'(E-W) by «50'(N-S). These data 
were obtained in the course of 12 ground-based observ- 
ing runs with various telescopes, as well as with the Wide 
Field Planetary Camera 2 (WFPC2) and the Advanced 
Camera for Surveys Wide Field Camera (ACS-WFC) on 
the Hubble Space Telescope (HST) . The accuracy of the 
zero-points of the optical catalog is believed to be better 
than 0.01 mag (see the bottom panel of Fig. 1). 

Initial photometry on individual images was performed 
with DA0PH0T IV, followed by simultaneous photometry 
for the 359 NIR and the 558 ground-based optical im- 
ages with ALLFRAME (Stetson, 1994). The HST images 
were reduced in a separate ALLFRAME run. The instru- 
mental IR magnitudes were transformed into the 2MASS 
photometric system. The standard error for one obser- 
vation of one star ranges from ~ 0.03 to ~ 0.06 mag, 
depending upon the filter (J,K) and the night. How- 
ever, with —22,000 individual observations of more than 
1,000 2MASS stars in J and -26,000 observations in K : 
we expect that our results should be on the 2MASS pho- 
tometric system to < 0.02 mag in the two bands. 

3. COMPARISON BETWEEN THEORY AND 
OBSERVATIONS 

The NIR CMDs present the obvious advantage of being 
minimally affected by uncertain and possibly differential 
reddening. However, image quality, angular resolution, 
and photometric precision have partially hampered the 
use of NIR photometry in crowded fields. Data plotted 
in the top panel of Fig. 2 show the K,J-K CMD based on 
data collected with NIR cameras on 4m telescopes (black 
dots, SOFIQNTT, ISPI@4m Blanco) and with MAD on 
VLT (red dots). The error bars plotted in the right side 
show that the photometric precision is of the order of a 
few hundredths of a magnitude over most of the mag- 
nitude range, increasing to — 0.1 mag at our adopted 
detection limit (K k, 20, J « 20.7). Note, in particular, 
that the MS in the low-mass regime shows a well defined 
knee. This feature is mainly caused by the collisionally 
induced absorption of H 2 at NIR wavelengths (Saumon 
et al. 1994). According to theory, the color and the shape 
of the bending depends on the metal content. However, 
the magnitude and color of the bend are, at fixed chem- 
ical composition, essentially independent of cluster age. 
This feature offers the unique opportunity to anchor clus- 
ter isochrones, and in turn to estimate the absolute age 
either as a color or as a magnitude difference between 
the bend and the cluster TO. The knee is a robust pre- 
diction, since for (M« 0.5-0.4 M Q ) the stellar structures 
are minimally affected by uncertainties in the treatment 
of convection, given that convective motions are nearly 
adiabatic (Saumon et al. 2008). In passing we note that a 
similar feature was already detected in two GCs (lu Cen, 
Pulone et al. 1998; M4, Pulone et al. 1999) and in the 
Galactic bulge (Zoccali et al. 2000) using deep and high 
angular resolution NIR (J,H) data collected with NIC- 
MOS@HST. To estimate the ridge line we represented 
each star detected in the J, if-bands by a Gaussian kernel 
with a sigma equal to the photometric uncertainty. The 
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individual Gaussians were summed and we identified the 
peaks along the MS. This ridge line was smoothed with 
a spline and the TO (red circle) and the MS knee (MSK, 
purple circle) were estimated as the points showing the 
minimum curvature in the two different color ranges (see 
also Table 1). 

To estimate the absolute age of NGC 3201, we have 
used the new set of isochrones for Y = 0.250, [Fe/H] 
= —1.50, and [a/Fe] ss 0.4 from the VandcnBerg et 
al. (2010, in preparation) compilation. [A constant en- 
hancement was not assumed for all a-elements; rather, 
the Cayrel et al. (2004) determinations of individual [el- 
ement /Fe] ratios in metal-deficient stars were adopted.] 
Fully consistent opacities and the latest nuclear reac- 
tions, along with treatments of H/He diffusion and tur- 
bulent mixing very close to those described by Proffitt 
& VandenBerg (1991) and Proffitt & Michaud (1991), 
respectively, were employed. Moreover, the free param- 
eters in their simple treatment of turbulence were set so 
that the predicted variation in the surface helium abun- 
dance during the main-sequence and subgiant phases 
agreed well with that reported by Korn et al. (2007). 
The latter had notable success explaining the observed 
abundances of a number of the heavy elements, using 
models that allowed for the settling and radiative ac- 
celerations of the metals; consequently, their prediction 
of the dependence of helium abundance on evolution- 
ary state represents the current best estimate. Although 
VandenBerg et al. have not treated the diffusion of the 
metals their isochrones are very similar to those used in 
the Korn et al. investigation, both morphologically and 
insofar as the predicted age at a fixed turnoff luminosity 
are concerned. 

The isochrones were transformed from the theoretical 
to the observed plane using the color-temperature (CT) 
relations derived by Casagrande, VandenBerg, & Stetson 
(2010, in preparation) from the latest MARCS model at- 
mospheres (Gustafsson et al. 2008). Small zero-point 
adjustments (0.03 mag) were applied to colors of the 
isochrones so that the synthetic colors for stars fainter 
than the turnoff agreed well with those derived from 
the IRFM-based, empirical CT relations for dwarf stars 
by Casagrande et al. (2009). Consistent fits of these 
isochrones to most of the CMDs for NGC 3201 that can 
be derived from BVRIJK S photometry indicated a clear 
preference for a true distance modulus ~ 13.35 ±0.11 
and a reddening E{B - V) « 0.24 ± 0.02. 

The extinction in the NIR bands was estimated using 
the McCall (2004) relations. The colored lines plotted in 
the top panel of Fig. 2 show the comparison with three 
isochrones at t =10,12 and 14 Gyr. However, the model 
isochrones in this CMD appear to be « 0.05 mag bluer 
than observed for magnitudes fainter than the MSK. 

To overcome subtle uncertainties in the cluster redden- 
ing, we performed a comparison between theory and ob- 
servations using a Wesenheit A-band 18 magnitude. This 
magnitude is reddening free, so the reddening uncer- 
tainty in the WK,J-K CMD occurs only in the horizon- 
tal direction. Data plotted in the bottom panel of Fig. 2 
show even more clearly the systematic drift in color of 
cluster isochrones when compared with observations. 

18 WK=K-A K /(Aj-A K ) [J-K] mag, where the Aj and A K arc 
the selective absorption according to McCall's relation. 



To further constrain the possible culprit (s) for this dis- 
crepancy we cross-correlated the NIR (J, A) photometry 
based on MAD images with the optical (V, I) photome- 
try based on ACS images. The top panels of Fig. 3 show 
the K,I - K (left) and the K,V - A (right) optical- 
NIR CMD. Data plotted in these panels display several 
interesting features, i)- The MS knee is still present, al- 
though less evident than in the NIR CMD. However, the 
approach we devised allows us to provide accurate mag- 
nitude and color estimates for the age indicators (sec 
Table 1). ii)- The same isochrones that in the NIR 
(K,J—K) CMD attain systematically bluer colors fainter 
than MSK provide — assuming identical values for dis- 
tance and reddening — a reasonable fit to both evolved 
and MS stars. However, the isochrones again attain sys- 
tematically bluer colors below the MS knee (A >18.5 
mag). The same outcomes apply if the comparison is 
performed using the reddening free Wesenheit A-band 
magnitude and the same optical-NIR colors (see the bot- 
tom panels in Fig. 3). 

As a further test of the isochrones we performed the 
same comparison for an optical I,V — I CMD. In par- 
ticular, we took advantage of the deep and precise pho- 
tometry from ACS and ground-based images. Data in 
Fig. 4 display a well defined cluster sequence ranging 
from the subgiant branch to the very low-mass regime 
(7 <~ 24, V - I ~ 2 mag). Note that a handful of can- 
didate cluster white dwarfs shows up in the left bottom 
corner (I ~ 24.5, V - I ~ 0.25 mag). The MS knee is 
less evident than in the optical-NIR CMDs, but both the 
high number of stars and the photometric precision per- 
mit a robust identification of the age indicators (see Ta- 
ble 1). Interestingly, we found that the same isochrones 
in the NIR and optical-NIR CMDs provide a good fit 
to the observations over the entire magnitude and color 
range. The isochrones are only slightly redder for mag- 
nitudes fainter than the MSTO. This clearly indicates 
that optical-NIR, and to a minor extent NIR, CMDs are 
an acid test of the precision of both evolution and atmo- 
sphere models. Even small changes in the input physics 
can be easily identified in the quoted CMDs. 

The comparison between theory and observations high- 
lights two relevant findings: i)— the adopted chemical 
composition is appropriate for NGC3201, since cluster 
isochrones simultaneously account for both optical and 
optical-NIR CMDs; ii)- the main reason in the discrep- 
ancy between theory and observations in the NIR and in 
the optical-NIR CMDs appears to be the A-band and to 
a small extent the /-band. The isochrones attain colors 
that are slightly either bluer (K,J-K) or redder (I,V-I) 
than the observed MSK. An independent and prelimi- 
nary test performed using the same isochrones, but a 
different set of CT relations (Casagrande et al. 2010) 
show a smaller color discrepancy and very similar cluster 
ages. 

4. RESULTS AND DISCUSSION 

The adopted cluster isochrones for ages ranging from 
10 to 14 Gyr indicate that the magnitude/age derivative 
of MSTO stars changes from AV/At- 0.093, to A7/At~ 
0.080 and to AA/A t ~ 0.062 mag/Gyr. This means 
that the age sensitivity decreases by «35% when mov- 
ing from the V to the A-band. On the other hand, the 
color/age derivative of MSTO stars is A(V - i")/At~ 
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0.013 for optical bands, A(V - JQ/At~ 0.032, A(I - 
K) /At- 0.019 for optical-NIR bands and A{J-K) /At- 
0.010 mag/Gyr for NIR bands. These color-age deriva- 
tives indicate that optical {V — I) and NIR ( J — K) colors 
attain similar sensitivities, while the optical-NIR (V — K, 
I — K) colors are, as expected, a factor of two-three 
more sensitive. A decrease of 0.2 dex in iron content 
implies a —20% decrease in magnitude/age derivatives, 
while an increase of the same amount causes minimal 
changes (<5%). The same changes in chemical composi- 
tion have minimal effects on the color/age derivatives. 

However, the reddening law by McCall (2004) gives 
selective absorptions of A//Ay-0.538 19 , A K /A V ~0.112 
mag and reddening corrections of E(V — K)/E(V — /)— 
1.923, E(I - K)/E(V - J)~ 0.923 and E(J - K)/E(V - 
J)~ 0.343 mag. These numbers indicate that the K mag- 
nitudes and J — K colors are respectively one order of 
magnitude and more than a factor of two less affected 
by reddening uncertainties (—10%) than the V — I colors 
and the I magnitudes. 

To constrain the error budget in the absolute age es- 
timates of GCs, based on the MSTO, we also need to 
account for uncertainties affecting the distance determi- 
nation. The uncertainty in the GC distance depends on 
the adopted standard candle and is still of the order of 
5% (Bono et al. 2008). Using the quoted error budget we 
estimated the absolute magnitude of the MSTO in the 
optical and in NIR bands (see Table 1). The error bud- 
get on absolute magnitudes is, as expected, dominated 
by the uncertainties in cluster distances. It is on average 
—0.11 mag for Mk, Mj and for the Wesenheit magni- 
tudes. The uncertainty of the chemical composition and 
the cluster isochrones will be discussed in a future paper 
(Dall'Ora et al. 2010, in preparation). Using the quoted 
isochrones we found absolute age estimates and uncer- 
tainties (see Table 1) for NGC 3201 that agree, within 
one sigma, with similar estimates in the literature. The 
weighted mean of the three age estimates based on the 
if -band is 11.48 ± 1.27, while that based on the WK- 
band is 11.55 ±1.53 Gyr. Optical and NIR (see Table 1) 
age estimates agree within one sigma. The relative ages 
(De Angcli et al. 2005) are typically affected by smaller 

19 To account for the difference in the effective wavelength be- 
tween the I-band adopted by McCall (A^^=7977 A) and by Lan- 

dolt (A^=8250 A) the selective absorption AijE{B — V) was 



errors, but they rely on predictions of evolved evolution- 
ary phases (HB, RGB) affected by different theoretical 
uncertainties. 

We also estimated the absolute ages of NGC 3201 from 
the difference in magnitude between the MSTO and the 
MSK using the same approach and the same isochrones. 
Values listed in Table 1 indicate that ages based on the 
new method agree, within the errors, with those based on 
the MSTO. The weighted mean of the three age estimates 
based on the if-band is 10.66 ± 0.58 and 11.63 ± 0.58 
Gyr based on the /-band, but the errors are more than 
a factor of two smaller. The same applies on average 
to the new age estimates based on the WK (10.82 ± 
0.68 Gyr) and on the VUT-band. The new age estimates, 
indeed, do not include uncertainties in the distance and 
reddening; their uncertainties are dominated by intrinsic 
photometric errors {<j,j : k —0.02, ay.i —0.01-0.02 mag). 

The new method to estimate the absolute age of 
GCs requires accurate NIR photometry «2-3 magnitude 
fainter than the TO region (Mj - M K «4-6). This 
means apparent magnitudes J — K «19-21 mag for 
a good fraction of GGCs and J — K «14-17 mag for 
Galactic open clusters. The former limiting magnitudes 
can be reached with 8-10m class telescopes, while the lat- 
ter are within the limits of 2-4m class telescopes. How- 
ever, the Extremely Large Telescopes and JWST will al- 
low us to apply this method to all the stellar systems 
in the Local Group — provided that suitable, local faint 
calibrating sequences exist. 
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changed, following Cardelli et al. (1989), from 1.714 to 1.652. 
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TABLE 1 

Magnitudes and colors of Main Sequence Turn-Off (MSTO) and Main Sequence Knee (MSK) in 

DIFFERENT OPTICAL AND NIR CMDs. 



CMD 


iriMSTO^ 


CImsto 




CI M SK h 


M M sto c 


M M sk c 


tMSTO d 


tMSTO-MSK C 


K, J-K 


16.26 


0.41 


18.30 


0.77 


2.83 


4.87 


11.8 ±2.3 


10.5 ± 1.0 


WK, J-K 


15.95 


0.41 


17.77 


0.78 


2.59 


4.41 


12.0 ±2.8 


10.4 ±1.1 


K, I-K 


16.26 


1.06 


18.31 


1.73 


2.83 


4.88 


11.8 ±2.3 


10.5 ± 1.0 


WK, I-K 


15.91 


1.06 


17.72 


1.63 


2.56 


4.37 


11.3 ±2.6 


10.9 ± 1.2 


K, V-K 


16.22 


1.96 


18.34 


3.10 


2.78 


4.90 


11.0 ±2.1 


11.0 ± 1.0 


WK, V-K 


15.92 


1.96 


17.81 


3.14 


2.57 


4.46 


11.4 ±2.6 


11.3 ± 1.2 


I, V-I 


17.28 


0.91 


20.46 


1.46 


3.53 


6.71 


11.5 ± 2.0 


11.6 ±0.6 


WI, V-I 


16.22 


0.91 


18.71 


1.43 


2.87 


5.36 


11.6 ±2.1 


11.2 ± 1.2 



a Apparent magnitudes of MSTO and MSK. 
b Color indices of MSTO and MSK. 
c Absolute magnitudes of MSTO and MSK. 

d Cluster age estimates (Gyrs) based on the absolute magnitude of MSTO. The uncertainties account for errors on 
photometry, reddening and distance. 

°Cluster age estimates (Gyrs) based on the difference in magnitude between MSTO and MSK. The uncertainties 
account for errors on photometry. 
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Fig. 1. — Top-Radial distribution of stars observed in the NIR (J,K) bands with different ground-based telescopes. Bottom-Same as 
the top, but for stars observed in the optical (V,I) bands with ground-based telescopes and with HST. 




Fig. 2.— Top - NIR (K, J - K) CMD of NGC3201. Black and red dots display stars observed with SOFIONTT, ISPIOBlanco and 
with MADOVLT. The colored solid lines display a set of cluster isochrones by vandenBerg et al. (2009), at fixed chemical composition, and 
different ages (see labeled values). These isochrones were transformed into the observational plane using the CT relations based primarily 
on MARCS model atmospheres (see the text). The adopted true distance and reddening are labeled. The large pink and purple circles 
mark the position of MS Turn-Off and of the MS knee. Bottom - Same as the top, but using the reddening free WK magnitude. 
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Fig. 3.— Top -K,I-K (left) and K, V — K (right) optical-NIR CMDs of NGC3201. Optical data were collected with ACS-WFC@HST, 
while the ii"-band data were collected with MAD@VLT. Symbols and lines are the same as in Fig. 2. Bottom — Same as the top, but using 
the reddening free Wk magnitude. 




Fig. 4.— Top - Optical I,V - I CMD of NGC3201. Data were collected with ACS-WFCOHST and ground-based telescopes. Symbols 
and lines are the same as in Fig. 2. The error bars on the right display intrinsic errors in magnitude and color. Bottom - Same as the top, 
but using the reddening free Wi magnitude. 



